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摘 要 : 沟渠 系统 “过 程 拦截 "是 现 阶 段 农业 非 点 源 污染 控制 和 管理 的 重要 手段 。 针 对 当前 农田 排水 沟渠 水 体 - 
底 泥 -植物 系统 内 各 介质 间 非 点 源 溶质 迁 移 转化 机 制 不 清 的 现状 , 在 简要 分 析 农 田 排 水 沟渠 水 - 底 泥 -植物 复合 
生态 结构 及 其 各 组 分 功能 特征 的 基础 上 ， 以 沟渠 中 水 体 为 中 心 ， 解析 了 非 点 源 溶 质 气 在 沟渠 系统 单一 介质 以 
及 不 同 介质 间 的 迁移 转化 过 程 ; 考虑 到 排水 沟渠 中 水 体 和 污染 物 均 以 沿 沟渠 方向 线性 迁移 为 主 应 用 水 流连 
续 方 程 和 非 保守 性 污染 物 迁 移 扩 散 方程 ,构建 了 农田 排水 及 其 中 非 点 源 污 染 物 在 排水 沟渠 中 的 一 维 迁移 模型 
并 结合 试验 沟渠 具体 情况 进行 了 适度 简化 ; 通过 内 业 可 控 性 监测 试验 ， 同 时 提出 了 沟渠 系统 水 - 底 泥 -植物 不 
同 介质 对 水 体 中 非 点 源 污 染 物 衰减 程度 影响 的 量化 计算 方法 ; 将 所 建 模型 与 计算 参数 应 用 于 野外 试验 沟渠， 
以 氨 态 所 和 硝酸 盐 所 为 例 进 行 模拟 ， 结 果 表 明 模 拟 结果 与 实测 值 结合 较 好 ， 其 中 氨 态 所 模拟 效率 系数 为 0.87， 
硝 态 氮 模 拟 效率 系数 0.93， 表 明 所 建 模型 及 模拟 手段 符合 试验 沟渠 实际 情况 ， 具 有 较 好 的 适用 性 。 
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Simulation of non-point Source Solute nitrogen transport in drainage 
ditches under water-sediment-aquatic plant synergy 


LI Qiangkun, HU Yawei, SU Xin 
(Institute of Yellow River Hydraulic Research, Yellow River Conservancy Commission, Zhengzhou 450003, China) 


Abstract: Along highlighted water environmental issues and gradual evolution of point source pollution (PSP) such as 
industrial waste, agricultural non-point Source pollution (AGNSP) due to non-scientific applications of fertilizers and 
pesticides has had an increasing attention in recent years. On the basis of Strengthening field “source control”, making the best 
use of agricultural drainage ditch/pond System (interception process) is critical for agriculture non-point source pollution 
control and management. Drainage ditches have compound ecosystem structures, which are consisted of water-sediment-plant 
synergy. Migration and transformation of non-point Source solute is important in the study of interception，control and 
management of agricultural non-point source pollution in drainage ditches. However, non-point source Solute transformation 
processes in each medium in water-sediment-plant System of farmland drainage ditch have remained largely unclear. Based on 


the analysis of ecological structure and drainage ditch characteristics, this paper described the transformation of non-point 
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source Soluble nitrogen in water-sediment-plant System of farmland drainage ditch. In combination with the characteristics of 


farmland drainage ditch water pollution migration, a model of agricultural water pollution migration in drainage ditch System 


was constructed. The model was based on continuity equation for flow and pollutants migration. It was used to test, monitor 


and advance a quantitative method of calculation of the effect of different water-sediment-plant media on attenuation degree of 


non-point source pollutant in water body. Using Qingshui drainage ditch in the People's Victory Canal Irrigation District in 


Henan Province as a case study, the model was used to calculate the parameters of a field experimental ditch loaded with 


ammonia nitrogen and nitrate nitrogen pollutants. The results showed that the Simulated and measured values fitted well. The 


simulation efficiency coefficient of ammonia nitrogen was 0.87 and that of nitrate nitrogen was 0.93. The results demonstrated 


that the model was available and effective for the simulation of drain pollution in ditches in agricultural farmlands. 


Keywords: Non-point source pollution; Farmland drainage ditch; Water-sediment-plant System; Nitrogen; Transportation 


simulation 


随 着 水 环境 问题 的 突出 以 及 工业 废 污 等 点 源 污 
染 的 逐步 整治 , 由 于 化 肥 农 药 不 科学 施用 而 引起 的 
农业 非 点 源 污染 问题 日 益 引 起 人 们 的 关注 和 重视 。 
结合 当前 我 国 农业 生产 现状 , 在 加 强 田 间 “ 源 头 控 
制 * 的 基础 上 ， 充 分 发 挥 沟渠 ( 塘 ) 系 统 “ 过 程 拦 截 * 是 
目前 农业 非 点 源 污染 控制 和 管理 的 重要 手段 ' “1。 灌 
区 网 状 分 布 的 大 量 排水 沟渠 为 非 点 源 溶质 净化 提供 
了 非常 有 利 的 “平台 ”。 农 田 狭 长 的 排水 沟渠 可 视 为 
县 有 河流 和 湿地 特征 的 独特 工程 化 生态 系统 中 ,在 
生态 学 和 物理 学 上 具有 与 线性 湿地 相似 的 功能 中。 
大 量 的 相关 试验 表明 农田 排水 沟渠 对 非 点 源 溶质 县 
有 较 好 的 拦截 效果 。Meuleman 等 5 指出， 天 然 沟 汇 
能 够 吸收 水 体 中 氮 磷 污染 物 ， 其 中 对 磷 素 的 去 除 率 
高 达 90%~95%; Kr6ger 等 "研究 发 现 , 一 条 长 约 
400 m， 宽 约 7 m 农田 排水 沟 对 可 溶性 无 机 氮 的 截留 
率 可 以 达到 57%， 对 磷 素 的 截留 率 可 以 达到 44%; 
吴军 等 外 在 江西 省 灌溉 试验 中 心 进行 的 沟 塘 湿地 对 
农田 面 源 污染 的 降解 试验 表明 : 沟 塘 湿地 对 农田 排 
水 中 总 氮 、 总 磷 的 综合 去 除 率 可 达到 70%~80%。 随 
着 沟渠 湿地 研究 的 持续 升温 ， 类似 结 果 还 有 很 
多 呈 ”总结 国 内 外 关于 农田 排水 沟渠 系统 净化 农业 
非 点 源 溶 质 的 相关 研究 ， 可 得 出 以 下 结论 : 四 农田 
排水 沟渠 对 非 点 源 溶质 县 有 明显 的 净化 效果 ; @ 应 
用 农田 排水 沟渠 “拦截 * 农 业 非 点 源 污染 尚未 形成 完 
善 的 理论 体系 ， 当 前 研究 主要 集中 于 局 部 或 单条 、 
多 条 沟渠 的 对 比试 验 、 水 生 植 物 备 选 及 生态 结构 形 
式 探讨 等 层面 ， 关 于 农田 排水 沟渠 水 - 底 泥 -植物 系 
统 内 各 介质 间 非 点 源 溶 质 的 迁移 转化 机 制 尚 不 清楚 ， 
关于 其 间 非 点 源 溶 质 的 运 移 模 拟 研 究 更 少 ， 这 也 是 


为 此 ， 本 项 研究 以 排水 沟渠 水 - 底 泥 -植物 协同 作用 
下 非 点 源 溶质 运 移 为 对 象 ， 基 于 水 流连 续 方程 和 污 
染 物 迁移 转化 方程 ， 以 监测 试验 为 基础 ， 通 过 估算 
水 - 底 泥 -植物 不 同 介 质 对 水 体 中 污染 物 衰减 程度 的 


影响 ,以 非 点 源 溶 质 氮 为 例 对 野外 试验 沟渠 中 氨 
态 所 和 硝酸 盐 气 的 运 移 进行 了 模拟 ， 以 供 后 期 进 一 
步 研究 借鉴 参考 。 
1 沟渠 生态 特征 与 营养 盐 迁 移 转化 
1.1 沟渠 生态 特征 

农田 排水 沟渠 系统 起 始 于 田间 毛 沟 或 农 沟 ， 经 
支 沟 、 干 沟 或 总 干 沟 排 入 外 界 大 面积 水 体 。 排 水 沟 
渠 的 生态 结构 一 般 可 分 为 3 个 层次 (图 1)。 上 层 为 水 
生 植 物 突出 水 面 之 外 的 荆 、 叶 ,并 随 排 水 沟渠 中 的 
水 位 高 低 有 所 变化 。 水 生 植 物 一 般 为 适 于 在 排水 沟 
渠 环 境 中 生长 的 萌 苇 (Phpragmaites communis)、 薄 草 
(Zpha orientalis) 等 挺 水 植物 ， 夏秋 季 生 长 旺盛 ， 冬 
春季 枯萎 ; 中 间 层 为 排水 沟渠 中 输送 的 农田 退 水 、 
浮游 生物 以 及 水 生 植 物 淹没 于 水 中 的 部 分 茎 、 叶 ， 
沟渠 中 水 位 高 低 随 农田 灌溉 (降雨 径流 ) 而 呈现 周期 
性 的 变化 ; 底层 为 水 生 植 物 的 根系 、 排 水 沟渠 中富 
合营 养 物 的 基质 底 泥 (淤积 物 ) 以 及 滋生 于 其 中 的 各 
类 微生物 。 


水 生 植 物 
Aquatic plant 
植物 根系 
Plant roots 

SA 
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水 体 Water body 


Deposited material 


1 农田 排水 沟渠 生态 结构 示意 图 
Schematic diagram of ecological Structure of 
farmland drainage ditch 


农田 排水 沟渠 水 - 底 泥 -植物 独特 的 生态 结构 决 
定 了 其 特有 的 生态 特征 。 排 水 沟渠 中 生长 着 适应 于 
此 环境 的 水 生 植物 ， 并 在 年 内 周期 性 地 生长 变化 ; 
渠 底 基质 淤积 物 随 水 位 升降 周期 性 地 暴露 、 济 没 ; 
底 泥 淤积 物 中 丰富 的 “营养 源 ” 保 证 了 水 生 植 物 的 生 
长 需求 和 其 中 各 类 微生物 的 繁衍 。 当 农田 排水 流 经 
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Fig. 1 
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时 ， 其 中 的 有 机 质 、 气 、 磷 等 营养 成 分 将 发 生 复 杂 
的 物理 、 化 学 和 生物 转化 作用 W"1。 
1.2 ”沟渠 气 营 养 盐 迁 移 转化 过 程 

由 于 过 量 施肥 灌 当 引起 的 田间 氮 素 流失 ,构成 
农田 排水 沟渠 中 主要 营养 盐 。 所 在 农田 排水 中 主要 
以 有 机 所 和 无 机 氮 两 种 形式 存在 ,农田 排水 沟渠 水 - 
底 泥 -植物 系统 中 气 的 迁移 转化 包括 矿 ( 氨 ) 化 作用 、 
硝化 与 反 硝 化 作用 、 植 物 吸 收 和 底 泥 吸 附 作 用 等 。 
具体 可 分 解 为 : 一 是 基质 底 泥 和 水 体 之 间 气 的 交换 ， 
包括 基质 底 泥 对 水 体 中 有 机 氮 的 沉积 吸附 以 及 基质 


底 泥 的 毛 态 氮 再 生 作 用 ， 基 质 底 泥 中 微生物 对 水 体 
中 氨 态 所 和 硝 态 氮 的 吸收 、 分 解 ,水 体 与 基质 底 泥 
中 氮 的 硝化 与 反 硝化 作用 等 。 二 是 水 生 植 物 和 水 体 
间 氮 的 转化 ,一 方面 水 生 植 物 可 以 直接 从 沟渠 水 体 
中 吸收 NH4、N0O; 离 子 并 同化 为 自身 所 需要 的 物质 
(蛋白 质 和 核酸 等 )， 从 而 净化 水 体 ; 另 一 方面 由 于 
植物 生长 的 周期 性 特征 ， 秋 冬季 节 植 物 枯 败 会 增加 
水 体 中 的 有 机 氮 ， 进 而 矿 化 为 氨 态 氮 ， 引 起 水 体 二 
次 污染 。 氮 在 农田 排水 沟渠 水 - 底 泥 -植物 系统 中 的 迁 
移 转化 可 概括 为 图 2。 
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2 农田 排水 沟渠 水 - 底 泥 -植物 中 氮 的 迁移 转化 


Fig.2 Migration of nitrogen in the system of water-sediment-plant of drainage ditch of farmland 


2 ”模型 构建 与 参数 确定 
2.1 沟渠 系统 水 体 及 污染 物 迁 移 模 型 

考虑 到 排水 沟渠 中 水 体 和 污染 物 均 以 治 沟渠 方 
向 纵向 迁移 为 主 基于 水 流连 续 方 程 和 非 保守 性 污 
染 物 迁 移 扩 散 方程 ， 文 献 [12] 将 排水 沟渠 中 水 体 及 
污染 物 运 移 描述 为 式 (1) 和 式 (2)。 


04 00 
一 十 一 二 1 
Ot Ox 4 ® 
OC 4 
人 人 oC 一 65 (2) 
ot 67 or Ot ot 


式 中 : 0 为 1 时 刻 通过 排水 沟渠 断面 A 的 流量 (ms ); 
4 为 侧 向 流量 强度 (考虑 沟渠 侧 渗 )m sm ); 4 为 
渠道 过 水 断面 面积 (m”); x 为 模拟 断面 距 控制 断面 的 
距离 (m); C 为 排水 沟渠 中 上 典型 污染 物 的 浓度 (mg:L ); 
Cs 为 各 排水 沟渠 “点 源 ” 汇 入 口 典 型 污染 物 的 浓度 
(mg '); C' 为 支 干 级 排水 沟 所 接纳 的 深层 排水 中 某 
ne ); 为 断面 平均 流速 (ms ); 
K; 为 污染 物 剪 切 流 弥 散 系数 (ms )， 又 称 纵向 弥散 
系数 ， 是 表征 流动 水 体 中 污染 物 在 泊 人 台 水 流 方 向 (或 纵 
向 ) 弥 散 的 速率 系数 ， 其 物理 意义 是 : 单位 时 间 内 污 


染 物 在 纵向 弥散 的 面积 ; 有 为 综合 衰减 系数 (4 ) ， 指 
污染 物 在 水 体 中 衰减 变化 的 速率 ， 其 物理 意义 是 单 
位 时 间 内 污染 物 在 水 体 中 衰减 掉 的 百分率 。 污 染 物 
综合 衰减 系数 不 仅 与 污染 物 可 降解 性 (或 易 衰 减 性 ) 
本 身 有 关 ， 还 与 水 动力 条 件 及 外 界 的 水 温 等 有 关 。 
本 项 研究 外 业 试验 沟渠 为 河南 省 人 民 胜 利 渠 
灌区 清水 渠 ， 试 验 段 长 度 600 m， 上 下 游 控制 断面 
闻 没 有 其 他 沟 ( 渠 ) 分 支 , 水 流 相 对 封闭 ， 上 下 游 变 化 
不 大 ; 同时 由 于 排水 沟渠 切割 较 浅 , 沟渠 深层 二 次 
排水 及 侧 渗 较 少 。 因 此 ， 本 项 研究 中 可 对 公式 (1)、 
公式 (2) 进 行 简化 ， 去 掉 侧 向 流量 强度 4g、 分 支 沟 ( 渠 ) 


汇 入 项 Cg、 沟 渠 深层 排水 项 C'， 进 而 简写 为 公式 
(3)、 公 式 (4)。 
A (3) 
Ot Ox 
oC 9C 62C 
(3 (4) 
式 中 : 流量 0、 断 面 面 积 4、 流 速 4、 典型 污染 物 浓 


度 C 参数 均 可 由 试验 测定 。 污 染 物 剪 切 流 弥 散 系 数 
KK, 可 以 利用 模拟 断面 上 的 实测 流速 或 浓度 资料 进行 
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推 求 (1。 
2.2 ”综合 衰减 系数 KK 的 推 求 

综合 衰减 系数 K, 不 仅 反 映 了 污染 物 在 水 体 作 
用 下 降解 速度 的 快慢 , 还 是 研究 区 域 水 质 污染 变 
化 、 计 算 水 环境 容量 以 及 纳 污 能 力 的 重要 参数 。 国 
内 外 关于 综合 衰减 系数 K 值 的 确定 已 有 较 长 的 研究 
历史 , 但 研究 方法 和 成 果 多 侧重 于 河流 和 湖泊 等 大 
面积 水 体 ， 对 农田 排水 沟渠 中 综合 衰减 系数 的 研究 
较 少 。 

结合 图 2 所 示 沟 渠 系统 内 部 氮 的 迁移 转化 过 程 ， 
进一步 研究 可 将 农田 排水 沟渠 衰减 系数 天 分 解 为 : 

K=K,tKntK, (5) 
式 中 : K, 指 水 体 自身 运动 引起 的 污染 物 豪 减 (d), 
如 水 体 表 层 蒸发 、 水 体内 部 发 生 的 硝化 与 反 硝 化 作 
用 等 ; Rn 为 沟渠 底层 淤积 物 引 起 的 水 体 中 污染 物 衰 
减 (4), 图 2 中 可 看 出 包括 有 机 和 氮 的 沉降 、 对 氨 态 氮 
的 吸附 和 再 生 等 ， 其 中 包括 底 泥 中 的 微生物 对 水 体 
中 污染 物 的 分 解 等 作用 ; K, 为 沟渠 中 水 生 植物 引起 
的 水 体 中 污染 物 衰减 (4 ,包括 植物 对 水 体 中 各 类 
合 氮 离子 的 吸附 净化 ( 正 值 ) 以 及 由 于 水 生 植 物 枯 败 
引起 的 水 体 中 有 机 和 氮 进 而 氨 态 氮 等 不 同 含 氮 污 染 物 
的 增加 ( 负 值 )。 

KKa 和 天 的 确定 可 采用 两 种 方法 ,一 是 必要 
的 监测 试验 ， 二 是 应 用 相关 动力 学 模型 进行 分 解 测 
算 。 本 项 研究 中 采用 监测 试验 的 方法 确定 。 

2.3 ”参数 确定 

为 确定 农田 排水 沟渠 中 K,、、K 和 有 ,不同 时间 
段 的 变化 ， 鉴 于 沟渠 外 业 试 验 多 种 影响 因素 的 不 可 
控 性 ， 选 择 在 黄河 水 利 科学 研究 院 节 水 与 农业 生态 
实验 基地 (河南 省 新 乡 市 ) 进 行 静态 监测 试验 。 

涉及 本 项 研究 的 试验 包括 3 个 相同 规格 的 PVC 
试验 模拟 箱 (长 x 宽 x 高 : 2 mx0.65 mx 0.65 m), 其 中 1 
号 模拟 箱 仅 加 入 水 样 ， 可 确定 K,; 2 号 模拟 箱 在 箱底 
均匀 铺设 20 cm 底 泥 , 试验 用 底 泥 取 自 项 目 外 业 实 
验 地 一 一 河南 省 人 民 胜 利 渠 灌区 清水 渠 ， 并 经 自然 
阴干 、 粉 碎 、 去 除 杂 质 ; 3 号 模拟 箱 在 铺设 20 cm 底 
泥 的 基础 上 ， 移 栽 已 经 培育 好 的 萌 苇 ， 初 始 种 植 密 
度 为 10 株 .m 。 试 验 时 间 段 为 2014 年 4 月 至 2015 
年 3 月 ， 试 验 过 程 中 结合 外 业 试 验 沟渠 实测 污染 物 
浓度 采用 水 体 中 添加 尿素 [CO(NH,);] 和 磷酸 氧 胺 
[INH4)2HPO4] 混 合 溶液 调控 水 体 中 所 、 磷 浓度 。 水 
样 每 10 d 取样 1 次 , 底 泥 样品 每 月 采集 1 次 。 各 项 
样品 提取 按照 相关 规定 进行 必要 的 前 处 理 后 进行 实 
验 室 化 验 分 析 ， 主 要 分 析 仪器 为 Smartchem 140 全 
自动 化 学 分 析 仪 。 


上 县 体 试验 布置 可 参见 文献 [14]。 
依据 水 样 、 底 泥 样 品 中 监测 数据 ,应 用 体积 法 
计算 水 体 中 含 氨 量 ， 见 式 (6)， 应 用 质量 法 计算 底 泥 
中 氮 含量 ， 见 式 (7)。 
W,=V,xC, (6) 
W, =Mx Co (7) 
式 中 : W,、 Wh 分别 为 水 体 中 、 底 泥 中 所 含 氮 素 (g); VV 
为 试验 箱 中 水 体 体 积 (m3)， 可 在 取样 时 进行 现场 测 
定 ; M 为 每 月 取样 时 试验 箱 中 底 泥 质量 (kg), 可 采 
用 装 箱底 泥 总 质量 扣除 每 次 取样 提取 的 底 泥 质量 获 
得 ; Cv 为 水 体 中 相应 氮 素 体积 比 浓 度 (mg:L 下 ); Cnm 为 
底 泥 中 质量 比 浓 度 (gkg )。 
依据 上 式 计算 结果 ,可 进一步 计算 K,, 见 式 (8)， 
Km、 可 同 理据 此 推算 。 


Ks=(W,i—Wyir,)/(W,ixa) (8) 
式 中 :i、i+l 指 相 邻 两 次 试验 取样 ; 4 为 两 次 取样 期 间 


间隔 天 数 。 

需要 指出 的 是 ,依据 上 述 各 式 计算 所 得 衰减 系 
数 是 应 用 静止 水 箱 模 拟 排水 沟渠 监测 试验 获得 , 与 
野外 农田 排水 沟渠 在 水 力 特性 、 曝 气 等 方面 还 有 一 
定 出 入 ， 其 优势 主要 在 于 实验 的 可 控 性 , 能 够 分 别 
估算 底 泥 、 水 生 植 物 对 水 体 中 溶质 的 影响 程度 。 
2.4 ”模型 输入 与 求解 

模型 计算 时 ,水 量 采用 第 1 类 (Dirichlet) 边 界 条 
件 , 排水 过 程 模 拟 上 边界 条 件 为 模拟 排水 沟 起 始 端 
水 量 ,下 边界 条 件 为 排水 沟 出 口 控制 水 位 ,计算 模型 
为 封闭 性 模型 ; 污染 物 计算 采用 第 2 类 (Neumann) 边 
界 条 件 ， 设 定 模 拟 排水 沟 起 始 端 污染 物 浓度 为 初始 
条 件 。 有 具体 求解 过 程 可 以 采用 非 耦合 求解 ， 即 先 单独 
求解 水 流连 续 方 程 , 求 出 有 关 水 力 要 素 后 ,再 求解 污 
染 物 迁 移 转 化 方程 ， 推 求 沟渠 典型 污染 物 输出 结果 。 

确定 模型 中 有 关 参 数 后 ,求解 时 水 流 方程 采用 
四 点 隐 式 差分 格式 离散 ,污染 物质 随 流 治 程 衰减 的 
剪 切 流 弥 散 方程 采用 显 式 差分 格式 ,应 用 计算 机 编 
程 进行 求解 。 
3 ”模型 应 用 与 效果 分 析 
3.1 ”模型 验证 

将 前 文 所 建 模型 及 相关 参数 应 用 于 项 目 外 业 试 
验 沟渠 河南 省 人 民 胜 利 渠 灌区 清水 渠 进行 验证 。 该 
渠 起 始 于 河南 省 新 乡 市 原 阳 县 ， 流 经 获 嘉 县 并 最 终 
进入 人 民 胜 利 渠 总 干渠 。 清 水 渠 因 早年 渠 水 清流 而 
得 名 , 近年 由 于 两 侧 农田 化 肥 、 农 药 施用 及 农田 排 
水 进入 而 使 水 质变 差 。 农 田 灌溉 期 间 ， 由 于 上 游 大 
量 排水 进入 ,两 侧 农 田 又 从 沟渠 中 提 水 进行 灌溉 ， 
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属于 灌 排 合 一 沟渠 。 试 验 沟 渠 段 布置 在 清水 党 与 东 模拟 过 程 中 ,以 图 3 中 治水 流 方向 1 号 断面 (D1) 
一 干渠 交叉 枢纽 下 游 约 500 m, 试验 段 总 长 度 600 m， ”所 测 氮 的 浓度 作为 输入 浓度 , 3 号 断面 (D3) 作 为 模拟 
相对 封闭 ,没有 其 他 进 水 口 和 排水 口 ， 治 水流 方向 ” 的 目标 浓度 ,以 氮 态 所 和 硝酸 盐 氮 为 例 进 行 模拟 。 
分 别 布设 1、2、3 号 断面 , 详 见 图 3 所 示 。 结果 见 图 4。 


野外 实验 区 示意 图 


Sketch of experimental area in the field 
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3 ”模型 验证 试验 布置 示意 图 


Fig.3 Schematic diagram of experiment layout for model verification 
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4 ” 氨 态 氮 (a) 和 硝 态 氮 (b) 实 测 数据 与 模拟 数据 (图 3 中 号 断面 的 输出 浓度 ) 对 比 图 
Fig.4 Contrast between simulated and measured (the output concentrations of the D3 section in the figure 3) concentrations of 
ammonia nitrogen (a) and nitrate nitrogen (b) 


3.2 ”模拟 结果 的 合理 性 分 析 实测 浓度 峰值 为 0.41 mg:L"'， 模拟 浓度 为 0.46 mg 

将 3 号 断面 实测 氨 态 氮 、 硝 态 氮 体积 比 浓度 分 别 。” ”实测 浓度 最 低 值 为 2014 年 8 月 12 日 的 0.10 mg 相 
点 绘 于 图 4。 可 以 看 出 图 4a 中 氮 态 氮 的 模拟 结果 与 ”应 模拟 值 为 0.07 mg-L”"， 峰 谷 值 略 有 出 入 ; 计算 模拟 时 
实测 值 在 浓度 峰值 和 低谷 有 差别 , 2014 年 4 月 24 日 ” 段 内 实测 氨 态 氮 体 积 比 浓度 平均 值 为 0.23 mg:L， 
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模拟 结果 平均 值 为 0.25 mg:L-， 二 者 较为 接近 。 比 
较 图 4a 中 实测 氨 态 氮 体 积 比 浓度 和 模拟 浓度 变化 曲 
线 ,二 者 变化 趋势 基本 一 致 。 图 4b 中 对 硝 态 氮 浓 度 
的 模拟 无 论 是 模拟 时 段 内 峰 谷 值 的 拟 合 程 度 以 及 浓 
度 曲 线 的 变化 趋势 ， 模 拟 结 果 均 好 于 图 4a 中 对 氨 态 
氮 的 模拟 , 计算 水 体 中 实测 硝 态 氮 体积 比 浓度 平均 
值 为 0.98 mg:L"',， 模拟 平均 值 为 1.02 mg:L"', 二 
更 为 接近 。 

为 进一步 分 析 模 拟 结果 的 合理 性 ， 同 时 采用 
Nash-Suttcliffe 模拟 效率 系数 (NSC) 对 模型 模拟 结 
进行 评价 。 计 算 公 式 见 式 (9): 

NSC=1 (Re 

(X= ZT) 

式 中 : Xos 为 实测 值 ， Xi 为 模型 计算 值 ，X ,为 实 
测 值 的 算术 平均 值 。 当 模型 计算 值 和 实际 监测 值 相 


(9) 


等 时 , NSC=1, 模拟 效果 最 好 ; 通常 NSC 在 0~1 之 间 ， 


NSC 越 大 ,说 明 计 算 值 与 观测 值 匹配 程度 越 好 。 

应 用 式 (9) 进 行 计算 , 氮 态 氮 模 拟 效 率 系 数 (NSC) 
为 0.87， 硝 态 氮 的 NSC 为 0.93， 同 上 述 分 析 结果 基 
本 吻合 ， 说 明 模 型 模拟 结果 均 可 以 接受 (大 于 0.75), 
其 中 模型 对 硝 态 氮 的 模拟 精度 略 高 于 对 和 氮 态 氮 的 模 
拟 精度 。 
4 讨论 与 结论 

目前 ,单纯 用 于 描述 农田 排水 沟渠 系统 水 体 及 
污染 物 运 移 的 数学 模型 较 少 ， 大 部 分 模型 都 是 在 湖 
泊 或 河流 模型 的 基础 上 发 展 起 来 的 。20 世纪 90 年 
代 ， 基 于 浅水 湖泊 富 营养 化 模型 PCLake 的 基础 上 ， 
答 兰 瓦 跟 宁 根 大 学 (Wageningen UR) 开 发 了 PCDitch 
模型 5…9， 该 模型 包涵 营养 物质 循环 、 东 类、 水生 
植物 、 底 栖 动 物 、 浮 游 动 物 和 鱼 类 等 模块 ， 同 时 包 
括 水 体 和 底 泥 区 域 ， 是 用 于 模拟 农田 沟渠 水 环境 变 
化 及 其 对 沟渠 生物 影响 的 富 营 养 化 模型 。 之 后 的 
Waterpas-models 模型 能 够 预测 沟渠 水 位 变化 和 耕 
作 方 式 等 对 农业 和 生产力、 农业 收益 和 环境 造成 的 影 
响 5， 模 型 模拟 发 现 ,过 高 的 沟渠 水 位 不 仅 会 减少 
农业 收成 和 降低 农民 的 经 济 收 益 ， 而 且 还 会 增加 地 
表 水 中 的 营养 负荷 ， 从 而 加 重 区 域 农业 非 点 源 污 
染 。 国 内 方面 ， 由 于 灌区 农田 排水 沟渠 系统 多 呈 
“ 树 ” 型 分 级 布设 ， 类 似 于 线性 河 网 系统 ,因此 沟渠 
系统 水 体 和 其 中 污染 物 的 模拟 多 借鉴 河道 水 质 水 
量 模型 进行 构建 ,水 量 模拟 采用 水 流连 续 方程 ， 其 
中 的 污染 物 采 用 一 维 的 非 保守 污染 物 迁 移 转 化 方程 
模拟 中 也 都 取得 了 不 错 的 模拟 效果 并 成 为 当前 


用 于 沟渠 系统 水 体 和 污染 物 运 移 模拟 的 主要 模型 。 

应 用 水 流连 续 方 程 和 非 保 守 污染 物 迁 移 转 化 方 
程 模拟 存在 的 最 大 问题 在 于 要 有 一 定 的 沟渠 试验 监 
测 资料 进行 明显 的 验证 和 校正 ， 其 中 的 难点 主要 在 
于 典型 污染 物 综合 衰减 系数 开 的 确定 。 由 图 2 中 排 
水 沟渠 水 - 底 泥 -植物 中 氮 的 迁移 转化 过 程 可 以 看 出 ， 
其 中 氮 的 迁移 转化 非常 复杂 , 文献 [18] 便 把 影响 综 
合 豪 减 系数 K 的 因素 分 为 6 类 ,并 认为 很 难 一 一 描 
述 , 采用 反问 题 求解 技术 确定 了 综合 衰减 系数 K。 文 
献 [19] 也 认为 综合 衰减 系数 K 是 农业 非 点 源 溶质 在 
排水 沟渠 系统 中 多 种 作用 的 综合 ， 其 中 包括 : 水 体 
中 污染 物 沉 降 、 已 经 沉降 污染 物 的 再 悬浮 、 渠 底 沉 
积 物 对 水 体 中 污染 物 的 截留 、 吸 附 等 多 个 过 程 ， 认 
为 在 实际 工作 中 这 些 过 程 、 作 用 很 难 监测 和 区 分 ， 
因此 一 般 取 其 综合 作用 效果 ,并 在 具体 计算 过 程 中 ， 
根据 不 同 级 别 排水 沟渠 典型 渠 段 断面 试验 资料 通过 
两 点 法 确定 了 综合 衰减 系数 K。 这 两 种 处 理 方法 实 
质 是 将 图 2 中 所 述 沟渠 系统 水 - 底 泥 -植物 复合 生态 
体系 中 各 类 迁移 转化 作用 进行 了 归 一 处 理 ， 不 能 分 
类 突出 各 生态 组 分 在 体系 中 的 作用 机 制 。 本 项 研究 
中 ,依据 排水 沟渠 水 - 底 泥 - 植 物 中 氮 的 迁移 转化 特 
征 ， 将 综合 衰减 系数 天 进一步 细 化 为 3 部 分 ， 分别 是 : 
水 体 自 身 运 动 引 起 的 污染 物 衰 减 系 数 K,、 沟 渠 底层 
淤积 物 引起 的 水 体 中 污染 物 豪 减 系 数 及 沟渠 中 
水 生 植物 引起 的 水 体 中 污染 物 衰减 系数 后， 并 结合 
模拟 实验 进行 参数 确定 ， 相 上 比 于 以 往 的 综合 作用 而 
言 研究 进一步 深入 ,也 更 加 符合 排水 沟渠 水 - 底 泥 - 
植物 单一 介质 及 不 同 介质 间 的 相互 作用 机 制 。 

从 模拟 结果 看 ， 文 献 [19] 应 用 式 (1)、 式 (2) 对 青 
铜 峡 灌区 第 一 排水 沟 沟 渠 系统 水 量 和 污染 物 进行 了 
模拟 ， 并 同样 采用 Nash-Suttcliffe 模拟 效率 系数 对 模 
型 模拟 结果 进行 了 评定 , 效率 系数 为 0.82; 本 次 研 
究 结 合 研 究 沟 渠 实 际 情况 ， 将 原 式 (1)、 式 (2) 简 化 为 
式 (3)、 式 (4)， 对 研究 沟渠 氨 态 氮 、 硝 酸 盐 氮 运 移 情 
况 进行 模拟 , 氨 态 氮 模 拟 效 率 系 数 为 0.87, 硝 态 氮 
模拟 效率 系数 为 0.93， 均 具有 较 好 的 模拟 效果 。 

综合 上 述 ， 本 文 研究 结论 可 归纳 如 下 : 

1) 在 简要 分 析 农 田 排水 沟渠 水 - 底 泥 - 植 物 独 特 
的 生态 结构 基础 上 , 提出 了 非 点 源 溶 质 氮 在 沟渠 水 - 
底 泥 -植物 系统 内 迁移 转化 过 程 。 

2) 结 合 农 田 排水 沟渠 水 污 运 移 特点 ， 基 于 水 流 
连续 方程 和 污染 物 迁 移 转 化 方程 ,初步 构建 了 农田 
排水 沟渠 水 污 运 移 模型 ， 并 通过 监测 试验 , 提出 了 
水 - 底 泥 -植物 不 同 介质 对 水 体 中 非 点 源 污 染 物 衰减 
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466 中 国生 态 农 业 学 报 2017 Chinaxiv 合 { 倪 鱼油 
程度 影响 的 量化 计算 方法 。 [10] 陈 海 生 , 王 光 华 , 宋 仿 根 ,等 . 生态 沟渠 对 农业 面 源 污染 物 
3) 将 所 建 模型 与 计算 参数 应 用 于 野外 试验 沟渠 ， 的 截留 效应 研究 [中 . 江西 农业 学 报 , 2010, 22(7): 121-124 
以 氨 态 气 和 硝酸 盐 氮 为 例 进行 模拟 ， 结果 表明 模拟 Chen H S，Wang G Song F G, et a Retention ng 
机 ee removal effects of ecological ditch on agricultural non-point 
结果 与 实测 值 结合 较 好 ， 其 中 氨 态 氮 模 拟 效率 系数 source pollutants[J]. Acta Agriculturae Jiangxi, 2010, 22(7): 
为 0.87, 硝 态 氮 模拟 效率 系数 0.93， 表 明 所 建 模型 121-124 
和 参数 计算 结果 县 有 可 活用 性 。 [11] 李强 坤 ， 胡 亚 伟 ， 孙 娟 . 农业 非 点 源 污染 物 在 排水 沟渠 中 
的 迁移 转化 研究 进展 [J. 中 国生 态 农业 学 报 , 2010,，18(1): 
参考 文献 References SO 
Li Q K, Hu Y W, Sun J. Migration and transformation of 
[1] 李 丽 华 , 李强 坤 . 农业 非 点 源 污染 研究 进展 和 趋势 [J]. 农 agricultural non-point source pollutants in drainage ditches[J]. 
业 资 源 与 环境 学 报 , 2014, 31(1): 13-22 Chinese Journal of Eco-Agriculture, 2010, 18(1): 210-214 
LiL H, Li Q K. The progress and trends of agricultural non- [12] 李强 坤 ， 胡 亚 伟 ， 孙 娟 ,等 . 基于 “ 源 ”"“ 汇 ”过 程 的 农业 非 
point source pollution research[J]. Journal of Agricultural 点 源 污染 模型 构建 及 应 用 [可 . 中 国生 态 农 业 学 报 ，2011， 
Resources and Environment, 2014, 31(1): 13-22 19(6): 1424-1430 
[2] 孙 宁 宁 . 水 稻 灌 区 新 型 农田 水 利 系统 防 污 减 污 试验 研究 [D]. Li Q K, Hu Y W, Sun J et al. Establishment and application 
武汉 : 武汉 大 学 , 2012: 80-95 of source/sink-based agricultural non-point source pollution 
Sun N N. Experiment and study of pollution mitigation by model[J]. Chinese Journal of Eco-Agriculture, 2011, 19(6): 
new irrigation and drainage System in rice-based irrigation 1424-1430 
districts[D]. Wuhan: Wuhan University, 2012: 80-—95 [13] 李强 坤 ， 李 怀 恩 . 农业 非 点 源 污染 数学 模型 及 控制 措施 研 
[3] Needelman B A, Kleinman PJ A, Strock J S, et al. Improved 究 : 以 青铜 峡 灌 区 为 例 [M]. 北京 : 中 国 环境 科学 出 版 社 ， 
management of agricultural drainage ditches for water quality 2010: 155-196 
protection: An overview[J]. Journal of Soil and Water Li Q K, Li HE.Mathematical Model and Control Measures of 
Conservation, 2007, 62(4): 171-178 Agricultural Non-Point Source Pollution[M]. Beijing: China 
4] Strock J S, Dell CJ, Schmidt J P. Managing natural processes Environmental Science Press, 2010: 155-196 
in drainage ditches for nonpoint source nitrogen control[J]. [14] 李强 坤 ， 宋 常 吉 ， 胡 亚 伟 ， 等 . 模拟 排水 沟渠 非 点 源 溶质 氮 
Journal of Soil and Water Conservation, 2007, 62(4): 188—196 迁移 实验 研究 [J]. 环境 科学 , 2016, 37(2): 520-526 
5] Meuleman A F M, Beltman B. The use of vegetated ditches Li Q K, Song C J, Hu Y W, et al. Transformation of non-point 
for water quality improvement[J]. Hydrobiologia, 1993, source soluble nitrogen in simulated drainage ditch[J]. 
253(1/3): 375 Environmental Science, 2016, 37(2): 520-526 
6] Kroger R, Holland M M, Moore M T, et al. Hydrological 15] Janse J H. A model of nutrient dynamics in shallow lakes in 
variability and agricultural drainage ditch inorganic nitrogen relation to multiple stable states[J]. Hydrobiologia, 1997, 
reduction capacity[J]. Journal of Environmental Quality, 2007, 342_343: 1_8 
36(6): 1646-1652 16] Janse J H. A model of ditch vegetation in relation to 
[7] Kriger R, Holland M M, Moore M T et al. Agricultural eutrophication[J]. Water Science and Technology, 1998, 37(3): 
drainage ditches mitigate phosphorus loads as a function of 139-_149 
hydrological variability[J]. Journal of Environmental Quality, 17] de Vos J A, van Bakel P J T, Hoving I E, et al. Waterpas- 
2008, 37(1): 107-113 model: A predictive tool for water management, agriculture, 
[8] 吴军 , 崔 远 来 , 赵 树 君 ,等 . 沟 塘 湿地 对 农田 面 源 污染 的 降 and environment[J]. Agricultural Water Management, 2006, 
解 试验 [J]. 水 电能 源 科学 , 2012, 30(10): 107-109 86(1/2): 187-195 
Wu J, Cui Y L, Zhao S J, et al. Experiment of agricultural [18] 王 超 , 王 沛 芳 , 唐 劲 松 , 等 . 河道 治 岸 芦苇 带 对 氨氮 的 削减 
non-point source pollution mitigation by field drain and pond 特性 研究 [可 . 水 科学 进展 , 2003, 14(3): 311-317 
wetland[J]. Water Resources and Power, 2012, 30(10): Wang C, Wang P F, Tang J S, et al. Ammonia sorption and 
107-109 retention characters in the riparian reed zone[J]. Advances in 
[9] 张 树 楠 ， 肖 润 林 ， 刘 锋 , 等 . 生态 沟渠 对 氮 、 磷 污染 物 的 拦 Water Science, 2003, 14(3): 311-317 
截 效 应 [J]. 环境 科学 , 2015, 36(12): 4516-4522 [19] 李强 坤 ， 胡 亚 伟 ， 李 怀 恩 . 农业 非 点 源 污染 物 在 排水 沟渠 


Zhang S N, Xiao R L, Liu F, et al. Interception effect of 
vegetated drainage ditch on nitrogen and phosphorus from 
drainage ditches[J]. Environmental Science, 2015, 36(12): 
4516-4522 


中 的 模拟 与 应 用 [J]. 环境 科学 , 2011, 32(5): 1273-1278 

Li Q K, Hu YY W, Li H E. Simulation and application of the 
agricultural non-point source pollutants in drainage ditch[J]. 
Environmental Science, 2011, 32(5): 1273—1278 


http:/www.ecoagri.ac.cn 


